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I. INTRODUCTION 
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We compare the Gibbs and Maxwell constructions for the hadron-quark phase 
transition in neutron and protoneutron stars, including interacting hyperons in the 
confined phase. We find that the hyperon populations are suppressed and that 
^Sj ■ neutrino trapping shifts the onset of the phase transition. The effects on the 

Q . (proto)neutron star maximum mass are explored. 
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(-H ' Some recent lattice QCD calculations have improved our knowledge on the nature of the 

CLc nucleon-hyperon (NY) and hyperon-hyperon (YY) interactions [Ij. Moreover, an intense 

Q ! experimental activity aimed at exploring hypernuclei has started at the J-PARC facility, 
which will hopefully clarify in the near future some unresolved questions on the hyperon 
interactions. This will help to improve the theoretical modeling of neutron star (NS) interi- 
ors, where hyperons are predicted to be present in large fractions, and many aspects of the 
equation of state (EOS) could be clarified. 
1> I There are essentially two microscopic approaches to derive the properties and the EOS 

of dense matter starting from the bare baryon-baryon interactions, namely the Brueckner- 
. Hartree-Fock (BHF) theory j2j and the variational method [3i] • Including the hyperon degrees 
of freedom, the EOS becomes very soft, and in the BHF approach the maximum mass for 
I a NS lies below currently observed values of about 1.7 Mq To overcome this problem, 
some possibilities have been suggested: i) The repulsive effects of three-body forces (TBF) 
among nucleons and hyperons, which may stiffen the EOS ii) Additional repulsion coming 
from hyperon-hyperon interactions (but no experimental data are available so far); iii) The 
hadron-quark phase transition jsl 0], which could lead to a partial suppression of the hyperon 



■ population, thus stiffening the EOS. This conclusion needs to be further explored, since at 



present there are many uncertainties regarding the quark matter EOS. 

The scope of this work is to present results on the composition and structure of neutron 
and protoneutron stars (PNS), which are produced in the aftermath of successful supernova 
(SN) explosions of very massive stars (M > 8 Mq). These objects can reach temperatures 
of the order of 30-40 MeV in their interiors, and are characterized by a temporary neutrino 
trapping with a conserved lepton fraction Yf, ^ 0.4 [sl, llf, lasting some seconds. Both thermal 
effects and neutrino trapping may result in observable consequences of the neutrino signa- 
ture from a supernova, and may also play an important role in determining whether a SN 
ultimately produces a cold NS or a black hole. In this paper, we focus on the hadron-quark 
phase transition, with the explicit inclusion of interacting hyperons in the confined phase. 
For hadronic matter we adopt the BHF microscopic EOS extended to finite temperature, 
and use the MIT bag model in the quark phase. 

This paper is organized as follows. In Section II we illustrate the Brueckner-Bethe- 
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Goldstone (BBG) many-body theory including hyperons at finite temperature. Sections III 
and IV briefly review the MIT bag model and the treatment of the hadron-quark mixed 
phase, respectively. In Section V we discuss our results regarding the structure of NSs and 
PNSs, in particular their maximum mass. Final conclusions are drawn in Section VI. 



II. HADRON PHASE: THE BBG THEORY AT FINITE TEMPERATURE 

In the recent years, the BBG perturbative theory has been extended in a fully micro- 
scopic way to the finite-temperature case , according to the formalism developed by Bloch 



and De Dominicis In this approach the essential ingredient is the two-body scattering 
matrix G, which in the finite-temperature extension is determined by solving numerically 
the Bethe- Goldstone equation, written in operatorial form as 

where the indices a, 6, c indicate pairs of baryons and the Pauli operator Q and energy E 
determine the propagation of intermediate baryon pairs. In a given baryon-baryon channel 
c = (12) one has 

g(i2) = [l-ni(A;i)][l-n2(A;2)], (2) 

-E(i2) = mi + m2 + ei(A;i) + e2(fc2) (3) 

with the single-particle energy ej(A;) = k^jlvtii + Ui{k), the Fermi distribution ni{k) = 
(^Q[ei{k)-iii]/T _j_ i)~^^ the starting energy W, and the two-body interaction (bare potential) V 
as fundamental input. The various single-particle potentials within the continuous choice 
are given by 

Ui{ki) = ReJ2 $^^(^2)(A:iA;2|G(i2)(i2)[^(i2)]|A:iA;2)^ , (4) 

2=n,p,A,S k2 

where ki generally denote momentum and spin. 

We choose the Argonne Vig nucleon-nucleon potential [l2| as two-body interaction V, 
supplemented by TBF among nucleons, in order to reproduce correctly the nuclear matter 
saturation point on 0.17 fm"^ E/A ?a -16 MeV. As TBF, we use the phenomenological 
Urbana model [l3|, which is actually reduced to a density- dependent two-body force by 
averaging over the position of the third particle (lol . 14 



Recently, the BHF approach has been extended by including consistently interacting 
hyperons at finite temperature |l5l| . In the hyperonic sector we employed the Nijmegen 
soft-core NY potentials NSC89 jl6|, fitted to the available experimental NY scattering data. 
Since at zero temperature only A and E~ hyperons appear in the neutron star matter up 
to very large densities [13], we restrict also the present study to these two hyperon species. 
Therefore, for fixed partial densities pi = ni{k), {i = n,p, A, S~) and temperature T, we 
solve self-consistently Eqs. ([1]) and (jlj) and calculate then the free energy density, which has 
the following simplified expression 



/ = E 



(5) 
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where 



(ni{k) \nni{k) + [1 - ni{k)] ln[l - ni{k)] 



(6) 



is the entropy density for component i treated as a free gas with spectrum ei{k). All 
thermodynamic quantities of interest can then be computed from the free energy density, 
Eq. ([5]); namely, the chemical potentials /Xj, pressure p, entropy density s, and internal energy 
density e read as 



P 



dl 
dpi ' 
2d{f/p) 



P 



dp 



df 



e = f + Ts 



(7) 
(8) 

(9) 
(10) 



where p = ^ • pi is the baryon number density. 

However, due to the large number of degrees of freedom (4 partial densities + tempera- 
ture), it is necessary to introduce some approximations in order to speed up the calculations. 
We adopt the so-called Frozen Correlations Approximation (FCA), assuming the correlations 
at T 7^ to be essentially the same as at T = 0. This means that the single-particle po- 
tential Ui{k) for the component i can be approximated by the one calculated at T = 0. 
Furthermore, we fit the numerical results by a sufficiently accurate analytical parametriza- 
tion. We find that the following functional form provides an excellent parametrization of 
the numerical data for the free energy density in the required ranges of nucleon density 
(0.1 fm"^ < pN < 0.8 fm"^), hyperon fractions (0 < pa/ptv < 0.9, < pe/ptv < 0.5), and 
temperature (0 MeV < T < 50 MeV): 



f{Pn, Pp, Pa, Ps, T) = FnPn 
+ (Fa + Faa + Fas) Pa ■ 

+ (Fe + FsE + Fsa) Pe 

with the parametrizations at zero temperature: 



C 



5/3 

-Pa 



2mA Ma 

C 5/3 



Fjv = (1 - 13) {aoPN + &oP??) + (3 {o-iPn + &ip5J) , 
Fy = (a^ + ayx + aYX^)pN + (&y + + h\x^)p'''' 



N 1 



Fyy' 
My 



O'YY'Pn Py' 
1+ [cy + Cyx) Pat 



(11) 



(12) 
(13) 
(14) 
(15) 



where pn = pn+Pp] x = pp/pN] /3 = (l-2a;)2; C = (3/5) (371^)2/3 ~ 5.742; and Y, Y' = A, S^. 
At finite temperature the expressions are extended as follows: 



= Fn{T = 0) 



+ 



aot Pn + {dot + eot ) In(pAf) + fot / pn 
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TABLE I: Fit parameters for the free energy density, Eqs. (jlltilSp . 



ao,bo,co,ai,bi,ci 


-286.6 397.2 


1.39 


88.1 


207.7 


2.50 


^A' '^A' ^A' "A' ^A' "A' '^A 


-403 688 


-943 


659 


-1273 1761 1.72 


n 1 QjOtI 7*? 


-114 





291 





1.63 


aAA,CAA,<iAA 


136 0.51 


0.93 








aAS,CAS,C?AS 














ass, CSS, c^ES 














asA,csA,rfsA 


89 0.33 


0.81 








^0 „1 „0 „1 
t-A, t-A, ^S' ^S 


0.22 -0.38 


-0.59 


-0.22 






ao, dp,eo, /o 


-202.0 396.9 


-190.6 


35.2 






ai,di,ei,fi 


-138.0 308.4 


-109.3 


31.2 






dA,e\, fA,gA,bA,c\ 


92.3 29.3 


39.4 152.3 


4.78 


3.95 


c^s,es,/s,5i;,fe,cs 


89.2 61.0 


63.6 


186.8 


1.13 


3.30 



+ ^dit^pN + (rfir + ei^) In(pw) + fitypN 
Fy = Fy{T = 0) 

+ {dyt^ + eyt^) In(pjv) + fytVpN + gyt^ Hpy 
My = My{T = {}) + hYt' ffN , 



(16) 

(17) 
(18) 



where t = T/(100MeV) and / and pi are given in MeVfm"^ and fm~^, respectively, (and 
mA,s in MeV"^fm"^). 

Technically, these parametrizations were obtained by performing about 10'^ BHF calcu- 
lations at zero temperature in the (p„, Pp, Pa, ps)^space and then using the FCA to generate 
finite-temperature results, increasing by about one order of magnitude the number of "data" 
points /(pn, Pp, Pa, Ps, ^)- The optimal values of the fit parameters, listed in Table I, were 
then determined hierarchically for cold nuclear matter, cold hypernuclear matter, hot nuclear 
matter, hot hypernuclear matter, so that the fits are optimized also in the more constrained 
cases 15 . 



III. QUARK PHASE: THE MIT BAG MODEL 



For the quark phase, we adopt the MIT bag model involving u, d and s quarks 0, 0, 18 
Probably, this model is too simple to describe quark matter in a realistic way, and we plan 
to adopt more sophisticated models in the future [l9-21|. We assume massless u and d 
quarks, s quarks with a current mass of = 150 MeV, and a density-dependent (but 
temperature-independent) bag constant. 



B{p) = + (5o - 5oo) exp 



^Po 



2i 



(19) 



with = 50 MeVfm-3, Bq = 400 MeV fm'^ and (3 = 0.17. 

This approach has been proposed in Ref. [gJ on the basis of experimental results on the 
formation of a quark-gluon plasma obtained at the CERN SPS. The above choice of the 
parameters allows the symmetric nuclear matter to be in the pure hadronic phase at low 
densities, and in the quark phase at large densities, while the transition density is taken as 
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a parameter. Several possible choices of the parameters have been explored in [6[, and all 
give a NS maximum mass in a relatively narrow interval, 1.4 Mc^ < M^ax ^ 1-7 Mq. For a 



more extensive discussion of this topic, the reader is referred to [l8j, where details regarding 



the MIT bag model at finite temperature are also given, which will not be repeated here. 



IV. MIXED PHASE 

With this parametrization of the density dependence of B we now consider the hadron- 
quark phase transition in (proto)neutron stars. We calculate the EOS of a conventional 
neutron star as composed of a chemically equilibrated and charge-neutral mixture of nucle- 
ons, hyperons, and leptons. We do not take into account anti-particles and muons in this 
paper, since their effects on the EOS are very small. 

We compare numerical results obtained using the Maxwell and the Gibbs constructions 
for the phase transition. We remind that in the Maxwell construction both phases are charge 
neutral and the conditions of mechanical and chemical equilibrium determine a coexistence 
phase, which is described by a constant pressure plateau typical of a liquid-vapor phase 
transition. On the other hand, in the Gibbs construction j22[ both the hadron and the quark 
phase are allowed to be separately charged, still preserving the total charge neutrality. This 
implies that neutron star matter can be treated as a two-component system, and therefore 
can be parametrized by two chemical potentials. Usually one chooses the pair (/ie,/i„), i.e., 
electron and baryon chemical potential. The pressure is the same in the two phases to 
ensure mechanical stability, while the chemical potentials of the different species are related 
to each other, satisfying chemical and beta stabihty. As a consequence, the pressure turns 
out to be a monotonically increasing function of the density, at variance with the Maxwell 
construction. We note that our Gibbs treatment is the zero surface tension limit in the 



calculations including finite-size effects |23l . |24 



The Gibbs conditions for chemical and mechanical equilibrium at finite temperature be- 
tween both phases read 

/in + /ie - ^J'Ve = f^d = fJ-s , (20) 
/ip + /ie - /ii., = /in = /iA = /i« + , (21) 

/is- + /ip = 2/i„ , (22) 
J9/^(/ie,/i„,T) = PQ(/ie,/X„,T) = • (23) 

From these equations one can calculate the equilibrium chemical potentials of the mixed 
phase corresponding to the intersection of the two surfaces representing the hadron and the 
quark phase, which allows one to calculate the charge densities and and therefore the 
volume fraction x occupied by quark matter in the mixed phase, i.e., 

XP? + (1 - x)pf = . (24) 

From this, the energy density bm and the baryon density pm of the mixed phase can be 
determined as 



Cm = X^Q + (1 - X)^H , 

pM = XPQ + (1 - X)PH ■ 



(25) 
(26) 
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V. RESULTS AND DISCUSSION 

In Fig. [1] we show the particle fractions obtained by applying the Gibbs construction to 
NS matter (T = MeV and x^^ = 0, upper panel), and to SN matter (T = 40 MeV and 
Ye = 0.4, lower panel). We observe that hyperons do not appear at all in the NS mixed 
phase, because the onset for the hadron-quark phase transition takes place at densities below 
the hyperon onset. In SN matter, thermal hyperons appear already at low density, but they 
are negligible. Hence, we conclude that the quarks suppress the hyperons in both cases. 

The resulting EOS p{pb) is displayed in Fig. [21 where the upper (lower) panel displays the 
calculations performed for NS (SN) matter. The pure hadronic (quark) phase is represented 
by the solid red (dashed green) curve, whereas the mixed phase is indicated by the blue 
(pink) broken line if the Gibbs (Maxwell) conditions are applied. We observe that the Gibbs 
mixed phase spans over a wide range of density, around 0.20-0.74 fm~^ for NS matter, 
and slightly less for SN matter, 0.29-0.71 fm~'^, because finite temperature counteracts the 



coexistence of phases, see Ref. [2J|. Compared to the EOS with the Gibbs construction, 
the density jump obtained with the Maxwell construction appears over a narrower density 
range in NS matter, whereas for SN matter the transition takes place at reduced pressure 
and ovei_a wider density range, which is mainly due to the trapping condition, as explained 
in Ref. 



The EOSs discussed above are the fundamental input for calculating the stable configu- 
rations of compact stars. For that, we use the well-known hydrostatic equilibrium equations 
of Tolman, Oppenheimer, and Volkov [25|. For completeness, as in [26], we have attached 
the Shen EOS [2^ in the low-density regime (p < 0.1 fm~^), since below this threshold 
clusterization sets in and nuclear matter becomes inhomogeneous. 

In Fig. Owe display the gravitational mass (in units of the solar mass, Mq = 1.99 x 10^^ g) 
as function of the central baryon number density Pb,c'- The upper (lower) panel shows 
results for NS (SN) matter, and the blue (pink) broken lines represent the calculations 
performed with the Gibbs (Maxwell) construction. For comparison, we also display the 
stable configurations in the pure hadronic phase (solid red line). 

We find that in NSs the pure hadronic EOS with hyperons gives rise to a maximum mass 
smaller than the canonical value, i.e. 1.44 Mq, whereas the hadron-quark phase transition 
allows to reach maximum masses slightly above 1.5 Mq (the difference between Gibbs and 
Maxwell case being very small). This increase is due to the fact that the hyperon populations 
are suppressed by the existence of quarks, as shown in Fig. [1] 

On the contrary, in the SN case all maximum masses are larger than 1.44 Mq, but here 
the phase transition lowers the maximum mass with respect to the purely hadronic matter, 
because trapping reduces the hyperon concentrations in the latter case compared to NS 
matter. Also, we observe that the value of the maximum mass with the Gibbs condition 
is slightly higher than the one with the Maxwell construction, which is due to the overall 
higher pressure of the matter in the former case, see Fig. [2J 

We also find that the stars with a mass close to the maximum one have quark cores. 



VI. CONCLUSIONS 



We have studied the quark-hadron mixed phase in cold neutron stars and hot neutrino- 
trapped protoneutron stars containing also hyperons, and compared explicitly the Maxwell 
and the Gibbs phase transition constructions. 
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We find that pure quark matter appears in the cores of compact stars in any situation and 
that the hyperon fractions are nearly completely suppressed by the appearance of quarks. 
Due to this reason the maximum NS mass is increased by the presence of quark matter. 
However, the simple MIT bag model used here is not capable to reach currently observed 
NS masses of about 1.7 Mq and it will be an important task for the future to implement 
more sophisticated quark models. 

Apart from this challenge, in this paper we did not consider finite-size effects in the mixed 
phase. In particular, it will be interesting to check the influence of trapped neutrinos on the 
pasta structures in future works. 
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FIG. 1: Particle fractions Xi = Pi/pB as a function of the baryon density pB for NS (upper panel) 
and SN matter (lower panel). The mixed phase is calculated by the Gibbs conditions. 
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FIG. 2: Pressure as a function of the baryon density for NS (upper panel) and SN matter (lower 
panel). The red (green) line displays the calculations for purely hadronic (quark) matter. The blue 
(pink) curve is the coexistence region when the Gibbs (Maxwell) construction is performed. 
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FIG. 3: Gravitational mass (in units of the solar mass) as a function of the central baryon density 
Pb,c for NS (upper panel) and SN (lower panel) matter. The purely hadronic configurations are 
represented by broken blue lines, whereas the red (green) lines represent the configurations of 
hybrid stars obtained applying the Gibbs (Maxwell) construction. 



